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CALCULATION OF THE DYNAMIC LONGITUDINAL STABILITY OF 
A TILT-WING V/STOL AIRCRAFT AND CORRELATION 
WITH MODEL FLIGHT TESTS 

By Joseph R. Chambers and Sue B. Grafton 
Langley Research Center 

SUMMARY 

An analytical study has been conducted to determine the important factors influ- 
encing the dynamic longitudinal stability characteristics of a typical tilt-wing V/STOL 
aircraft. Calculations have been made for the initial condition of steady level flight at 
wing incidences corresponding to speeds ranging from hovering to conventional forward 
flight. The results of the calculations have been compared with qualitative measure- 
ments of dynamic stability obtained during free-flight tests of a 1/9-scale model of the 
aircraft. 

The results of the investigation indicated that the control-fixed motions of the air- 
craft without artificial stabilization in hovering flight would be dominated by an unstable 
oscillation similar to that displayed by most helicopters. As the transition to conven- 
tional forward flight progressed, stability characteristics were encountered in which 
aperiodic divergent modes of motion, as well as unstable oscillations, were present. The 
conventional-aircraft short-period and phugoid oscillations began to appear at the high- 
speed end of the transition. In general, the analytical results agreed with the motions 
observed during the free-flight model tests. These results also indicated that the unstable 
oscillation occurring in the hovering and low-speed flight regions can be stabilized by the 
addition of a combination of pitch-rate and pitch -attitude stabilization, but that angle-of- 
attack stability must be increased if the aperiodic divergent modes of motion are to be 
made stable. 


INTRODUCTION 

A considerable amount of data relating to qualitative measurements of the dynamic 
stability characteristics of dynamically scaled models of tilt -wing V/STOL aircraft has 
been published by the National Aeronautics and Space Administration. (For examples, 
see refs. 1 to 3.) At the present time, however, there is a lack of detailed information on 
the modes of motion and on the effect of various stability derivatives on the stability char- 
acteristics of this type of aircraft. The deficiency has been due in large part to the 



unavailability of representative values of the static and dynamic stability derivatives 
for V/STOL aircraft. 

As the first step of a program to provide some basic detailed information regarding 
the dynamic stability characteristics of tilt -wing configurations, static and dynamic lon- 
gitudinal stability derivatives were determined for a powered 1/9-scale model of a four- 
propeller tilt-wing V/STOL transport (ref. 4). Prior to the force tests, the model was 
flown in free -flight tests during which qualitative measurements were made of the 
control-fixed longitudinal motions at several speeds in the transition range (ref. 3). The 
experimentally determined values of the stability derivatives were the necessary inputs 
for an analytical investigation of dynamic stability, and the model flight tests provided 
data for purposes of correlation and validation. The present investigation was undertaken 
(1) to determine appropriate methods with which to analyze the longitudinal dynamics of 
tilt -wing vehicles, (2) to calculate the dynamic stability of a typical tilt -wing aircraft 
over the transition flight range, (3) to correlate the results of the calculations with those 
of the model free-flight tests, and (4) to determine the effects of the various static and 
dynamic stability derivatives on control-fixed dynamic stability. 

SYMBOLS 

All stability derivatives are presented with respect to the body system of axes 
shown in figure 1. Inasmuch as conventional nondimensional coefficients lose their sig- 
nificance and tend to become infinite as airspeed approaches zero, the stability deriva- 
tives are presented in dimensional form. The units used for the physical quantities in 

\ 

this paper are given both in the U.S. Customary Units and in the International System of 
Units (SI). Factors relating the two systems of units are given in reference 5. 

A,B,C,D,E coefficients defined in appendix A 

mean aerodynamic chord, ft (m) 

cycles required for oscillation to damp to one-half amplitude 
cycles required for oscillation to double amplitude 
force along X body axis, lb (N) 
force along Z body axis, lb (N) 
acceleration due to gravity, ft/sec^ (m/sec2) 
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tail incidence, deg 
wing incidence, deg 

moment of inertia about Y body axis, slug-ft 2 (kg-m2) 

aircraft linear dimension, ft (m) 
model linear dimension, ft (m) 

value of lift for longitudinal acceleration equal to zero at an angle of attack 
of 0°, lb (N) 

mass, slugs (kg) 

pitching moment, ft-lb (m-N) 

period of oscillation, sec 

pitching velocity, rad/sec 

dynamic pressure, pV 2 /2, lb/ft 2 (N/m2) 

Laplace operator, a + jco, 1/sec 

wing area, ft 2 (m2) 

time required for a mode of motion to damp to one-half amplitude, sec 
time required for a mode of motion to double amplitude, sec 
perturbation velocities along X and Z body axes, ft/sec (m/sec) 
trim velocity, ft/sec (m/sec) or knots 

ratio of amplitude of perturbation velocity along Z body axis to perturbation 
velocity along X body axis 
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V 

AV 

w/s 

X,Y,Z 

a 

Aar 
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ratio of amplitude of perturbation in pitch angle to perturbation velocity 

deg / deg \ 


along X body axis. 


ft/sec Im/seci 




velocity, ft/sec (m/sec) 

velocity increment, ft/sec (m/sec) 

wing loading, lb/ft^ (N/m2) 

body reference axes (fig. 1) 

angle of attack, deg or rad 

increment of angle of attack, rad 

flight-path angle, positive for climb, rad 

ratio of damping present in oscillatory mode of motion to value required for 
critical damping 

pitch angle, positive when nose is above horizon, deg or rad 

scale factor, l /l M 

air density, slugs/ft^ (kg/m3) 

real part of root of characteristic equation, 1/sec 

imaginary part of root of characteristic equation, rad/sec 


OO 


n 


undamped natural frequency of oscillatory mode, rad/sec 
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Dimensional stability derivatives: 
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Dot over a symbol indicates a derivative with respect to time. 

METHOD OF ANALYSIS 

The initial phase of the study consisted of an analytical determination of the dynamic 
longitudinal stability characteristics of a typical tilt -wing V/STOL aircraft for level flight 
(values of a and y initially equal to zero) at several values of wing incidence. A large 
part of the investigation was devoted to arriving at suitable methods of analysis — in par- 
ticular, to a determination of whether the classical linearized, small-perturbation equa- 
tions of motion normally used in conventional aircraft stability analysis would adequately 
describe the control-fixed motions of a tilt-wing vehicle. The validity of linear equations 
of motion when applied to tilt -wing aircraft might be questionable because (1) large- 
amplitude control-fixed motions are expected to occur for V/STOL aircraft (see refs. 1 
to 3) and (2) aerodynamic nonlinearities have been reported for tilt-wing configurations 
over part of the transition flight range (see ref. 4). The classical approach is, of course, 
the simplest method of analyzing dynamic-stability problems. The results of the calcula- 
tions made with the linear equations of motion were compared with the scaled-up results 
of the free-flight model tests to determine the validity of the calculations. In conjunction 
with the linear analysis, values of the various longitudinal stability derivatives appearing 
in the equations of motion were varied not only to determine which derivatives had signifi- 
cant effects on dynamic stability and must therefore be known with a reasonable degree of 
accuracy, but also to determine which derivatives were so insignificant that they could be 
neglected for purposes of analysis. 
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Description of Aircraft 


The aircraft used for the calculations was a four-propeller, tilt-wing V/STOL 
transport and is believed to be typical of present-day flapped tilt -wing designs. A three- 
view sketch showing the more important dimensions is presented in figure 2. The con- 
figuration employed counterrotating propellers, which resulted in cancellation of the pro- 
pulsive gyroscopic effects and, in turn, in uncoupling of the longitudinal and lateral equa- 
tions of motion. The mass and geometric characteristics of the V/STOL transport are 
given in table I. The values were obtained by scaling up the corresponding values for the 
1/9-scale model as measured during the free-f light tests of reference 3. Scaling factors 
used are given in table n. The horizontal and vertical locations of the center of gravity 
were assumed to vary with wing incidence angle in the same manner as those for the 
free-flight model; the variations are presented in figure 3. The configuration employed 
a full-span flap and an all-movable horizontal tail. These two surfaces were programed 
to deflect with wing incidence. The programed variations of the flap deflection and 
horizontal-tail incidence (identical to those of the free-flight model) are shown in figure 4 
and are the same as those used in reference 4. Unbalanced pitching moments occurring 
with this particular programed variation were assumed to be trimmed by an auxiliary 
device (the free-flight model was trimmed with a jet-reaction type of control) with no 
resulting effect on the longitudinal stability derivatives. 

The results of reference 4 indicate that the tail-incidence program used produced 
excessive values of tail incidence at the higher transition speeds — that is, large nose- 
down pitching moments were produced by the horizontal tail. This point should be kept 
in mind inasmuch as the tail-incidence program produced negative (statically destabi- 
lizing) values of the velocity stability derivative M u . In other words, when the speed of 
the model was increased from the trim flight condition, the pitching-moment contribution 
of the jet-reaction-control trimmer remained unchanged, but the moment increment due 
to the horizontal tail became more negative. The result is that different trimming pro- 
cedures for tilt-wing configurations may lead to completely different values of the longi- 
tudinal stability derivatives. For example, the velocity stability derivative M u will 
vary at higher speeds, depending upon whether the vehicle is trimmed with the horizontal 
tail or another trimming device such as a tail rotor. 

Correlation of the results of the present study with those of similar configurations 
(refs. 6 and 7, for example) was not attempted because of differences in tail-incidence 
programs and mass and inertial characteristics. 

Stability Derivatives 

The static -force -test data of reference 4 (obtained for conditions of constant power) 
were linearized by measuring representative slopes of the force and moment data plotted 
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as functions of angle of attack and velocity. Slopes were obtained for forward flight for 
each wing incidence at the power condition that gave zero net longitudinal force at an 
angle of attack of zero (Fx = 0 at a = 0°). The static stability derivatives were then 
expressed in terms of the perturbation velocities u and w by the approximate rela- 
tions AV ~u and U 0 Ao! ®w. The change of variables was desirable inasmuch as the 
values of the stability derivatives for hovering (where a is undefined) and for forward 
flight then become directly comparable. This procedure was not required for estimation 
of the derivatives for hovering flight inasmuch as the data presented in reference 4 were 
measured as direct functions of u. Values of the stability derivatives representing rate 
of change of vertical and horizontal forces with respect to the variables in the equations 
of motion were divided by the value of the model mass, whereas derivatives of pitching 
moment were divided by the model moment of inertia in pitch. The dimensional stability 
derivatives were then scaled up to full-scale values by multiplying the model values by 
the appropriate scaling factors given in table II. Multiplication of model values by the 
scaling factors results in no change in dynamic stability — that is, the model and aircraft 
have the same values for the nondimensional stability characteristics ^for example, 

Ci /2 and 

Several difficulties arose in the estimation of the stability derivatives. First, the 
results of reference 4 did not report data for the dynamic derivatives in hovering flight. 
Preliminary calculations indicated that the only dynamic derivative which could signifi- 
cantly affect the hovering dynamics was the damping-in-pitch derivative Mq. Unpub- 
lished forced oscillation tests subsequent to those of reference 4 indicated that the model 
had values of Mq equal to zero in hovering flight. The small value of the derivative 
was also confirmed by the results of damping-in-pitch tests made with a similar config- 
uration, as reported in reference 6. For the foregoing reasons, the dynamic derivatives 
for hovering flight were assumed to have values of zero for the aircraft without artificial 
stabilization. A second difficulty which arose during the estimation of stability deriva- 
tives was that the vertical damping derivative Z w was not measured during the model 
force tests. For hovering flight, the value of Z w was also assumed to be equal to zero. 
This assumption was not a critical one for the dynamic -stability calculations because the 
vertical translatory mode of motion is normally uncoupled from an unstable oscillation 
involving horizontal translation and pitching motion which dominates the control-fixed 
motions in hovering flight (see refs. 1 to 3). Information regarding height control of 
hovering VTOL aircraft (where the value of Z w may be relatively important) can be 
found in references 8 and 9. A third difficulty arose when the static -force -test data of 
reference 4 showed significant nonlinear variations of pitching moment with angle of 
attack for a wing incidence of 25°. (The static -force -test data for this condition are pre- 
sented in fig. 5 for reference.) For the trim condition ^L 0 /qooS = 7.0 ), the model was 

statically stable with respect to angle of attack (negative values of Mq) for small negative 
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angles of attack but was unstable (positive values of Mo) for small positive angles of 
attack. It might be expected that this nonlinearity would preclude the valid usage of linear 
equations for that particular wing incidence. In order to evaluate the validity of the linear 
equations, slopes were measured at both small values of positive and negative angles of 
attack. Two values of angle -of -attack stability were therefore obtained for i w = 25°. A 
final difficulty was that the nature of the forced oscillation tests of reference 4 led to 
dynamic -stability parameters that were combinations of stability derivatives due to pitch 
rate (q) and rate of change of angle of attack (<i). For purposes of the calculations, the 
measured parameters were assumed to be due entirely to pitch rate (for example, the 
measured parameter Mq + Mq» was used as Mq). 

The stability derivatives as determined by the foregoing methods and used in the 
calculations are presented in figure 6 and in table HI. 

Method and Scope of Calculations 

Calculations were made to determine the longitudinal dynamic stability character- 
istics of the vehicle for level flight at wing incidence angles of 90° (hovering), 65°, 50°, 
25°, and 10°. The calculations were not intended to describe the motions of an accel- 
erated transition in which the stability derivatives are functions of time. Instead, the 
investigation was made to determine the stability of the vehicle in steady flight at each 
value of wing incidence. Presented in figure 7 and table m are the full-scale trim veloc- 
ities for the particular wing-incidence values. The linearized, small-perturbation equa- 
tions of motion used for the calculations are given in appendix A. Further simplification 
of the equations for hovering flight is discussed in appendix B. The stability parameters 
obtained from the solutions of the linear equations were in the form of damping and peri- 
odic characteristics of the modes of motion. The damping characteristics included time 
to damp to one-half amplitude (ti/2)> damping ratio (£), and number of cycles required to 
damp to one-half amplitude (Ci/2)- Positive values of these parameters indicate sta- 
bility and negative values indicate dynamically unstable modes of motion. The periodic 
characteristics were the damped period (P) and the undamped natural frequency (o> n ) of 
oscillatory modes. Additional calculations were made to determine the general nature of 

the modes of motion as expressed by the amplitude ratio j^-j. Small values of this 

parameter indicate modes of motion containing little change in angle of attack (for example, 
the conventional-aircraft phugoid oscillation). Large values indicate large variations in 
angle of attack (similar to the conventional-aircraft short-period oscillation). 

In conjunction with the linear analysis, the root-locus technique of reference 10 
was used to illustrate graphically the variation of the roots of the characteristic equation 
with changes in the various stability derivatives. This method presents the path, or locus, 
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of the solutions on the complex plane as the value of an individual stability derivative is 
varied. (Additional information on the application of the method can be found in ref. 11.) 
Presented in figure 8 are the features of the complex plane as applied to dynamic sys- 
tems. For stability, all roots of the characteristic equation must have negative real 
parts (a negative). This requirement means that for stability all roots must be located 
in the left half of the complex plane shown in figure 8. The imaginary axis (at a = 0) is 
therefore a line of neutral dynamic stability for both real and complex roots. Lines 
representing constant values of time to halve (or double) amplitude are parallel to the 
imaginary axis (ct = Constant). Lines representing constant values of damped period are 
parallel to the real axis (ju> = Constant). Radial lines emanating from the origin repre- 
sent constant values of damping ratio or cycles to halve (or double) amplitude 
(w/ct = Constant). Circles with centers at the origin represent lines of constant undamped 

natural frequency (\Ja^ + co^ = Constant) . 

Although calculations were made for the nonlinear condition of i w = 25°, the 
results were rather meaningless for the present configuration inasmuch as the motion of 
the vehicle depended on the direction of the disturbance. Additional calculations were 
therefore made by using the basic force-test data of figure 5 in a digital computer pro- 
gram. This program utilized complete nonlinear equations of motion which were solved 
by numerical methods to produce representative time histories of motions obtained at 
i w = 25°. 


RESULTS AND DISCUSSION 
Calculated Stability of Aircraft 

Hovering flight .- The roots describing the modes of motion on the complex plane 
are plotted in figure 9 for hovering flight. A real root with a value equal to zero (a = 0) 
representing the neutral stability of the vertical motion of the vehicle appears at the 
origin. A negative (stable) real root associated with an aperiodic mode involving both 
horizontal translation and pitching motion is shown on the negative real axis. The addi- 
tional roots were a complex pair whose real part a was positive (unstable). As pointed 
out in appendix B, the oscillation involved both fore-and-aft motions and pitch-attitude 
changes. After a disturbance from hovering flight, the damped aperiodic mode would 
rapidly subside and the dynamic stability of the vehicle would be dominated by the unstable 
oscillation. 

The calculated stability characteristics of the aperiodic and oscillatory modes cor- 
responding to the roots are presented in table IV. The vertical translatory mode is 
neutrally stable (tj /2 = °°), because the derivative Z w was assumed to have a value of 
zero, as previously discussed. The other aperiodic mode is relatively highly damped 
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1/2 = 0.8 secj. The oscillation is unstable in terms of time required to double ampli- 
tude (t2 = 2.4 sec) with a moderately long period (P = 9.7 sec). These modes of motion 
are similar to those displayed by most helicopters in hovering flight, although the oscil- 
lation is more unstable for the tilt -wing aircraft. For example, the basic helicopter used 
in the investigation of reference 12 also possessed an aperiodic mode (ti/2 = 0.31 sec) 

and an unstable oscillation (P = 8.9 sec, t2 = 6.6 sec) in hovering flight. The point to be 
inferred is that dynamic instability can be tolerated for tilt -wing aircraft in hovering 
flight if they are operated in the same manner as helicopters. The fact that complete 
dynamic stability is not required for satisfactory handling qualities is also demonstrated 
by conventional aircraft with unstable phugoid or spiral modes of motion. 

Transition flight.- The data of table IV show that as the wing incidence is reduced 
from 90° to 65°, the most noticeable change from the stability characteristics of hovering 
flight is the appearance of an unstable aperiodic mode of motion. Analysis revealed that 
this instability of the vehicle was a result of angle -of -attack instability ^positive values 

of M w ) as indicated by the data of figure 6(a). The data of table IV also show that as 
the wing incidence was reduced from 90° to 65°, the unstable oscillation increased in 
period (from P = 9.68 sec to P = 14.11 sec) and was less unstable in terms of time 
required to double amplitude (t2 = 4.29 sec). The stable aperiodic mode for i w = 65° 
had about the same value of tj /2 as f° r hovering flight. The data of table IV also indi- 
cate that the aperiodic modes of motion involve larger angle-of-attack changes than 
does the oscillatory mode. 

As the transition progresses and the wing incidence is reduced to 50°, the larger 
value of angle-of-attack instability (see fig. 6(a)) results in greater instability of the 
unstable aperiodic mode. The oscillation is still unstable and continues to increase both 
in period (P = 19.40 sec) and in time required to double amplitude (t2 = 16.50 sec). The 
general nature of the aperiodic and oscillatory modes for i w = 50° is apparent in that 
the aperiodic modes contain larger changes in angle of attack and pitch angle than does 
the oscillatory mode. 

Further reduction in wing incidence to 25° results in the aerodynamic nonlinearities 
previously discussed. The results of calculations from the linearized values of the static 
data based on positive and negative angle-of-attack disturbances indicate that the motion 
initiated by positive angle-of-attack disturbances would be dominated by an aperiodic 
divergence (primarily in angle of attack and pitch angle) caused by the static instability 
with respect to angle of attack. For negative angle-of-attack disturbances, on the other 
hand, two oscillatory modes similar to the classical phugoid and short-period oscillations 
are present. The results of the linearized calculations for negative angle-of-attack dis- 
turbances are meaningless in that, as the aircraft angle of attack returns to the trim 
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condition and overshoots, the static instability will prevail — that is, the aircraft will 
always diverge in a nose-up direction. The results of additional calculations made by 
using automatic -computing equipment and the basic nonlinear data of figure 5 are pre- 
sented in figure 10. These data are time histories of the motions resulting from initial 
angle -of -attack disturbances of 5° and -5°. The motions for positive angle-of-attack 
disturbances show simply that the model diverges. For negative angle-of-attack dis- 
turbances, the vehicle initially tends to return to the trim condition; but, as the angle of 
attack overshoots, the model diverges as it did for the positive disturbance. The results 
for i w = 25° indicate that the classical linearized equations of motion are not applicable 
to some of the normal flight conditions of tilt-wing aircraft. 

For the lowest value of wing incidence of this investigation (i w = 10°), the results 
show a highly damped relatively short period oscillation, a stable aperiodic mode, and an 
unstable aperiodic mode. The amplitude ratios presented in table IV for i w = 10° indi- 
cate that the oscillation primarily involved changes in angle of attack and pitch angle 

^large values of j^-J and whereas the aperiodic modes primarily involved large 

changes in forward speed ^relatively small values of j^-j and |-^jj . The instability of 

one of the aperiodic modes was found to be a result of static instability of the aircraft. 

In this connection, it should be pointed out that the static stability of a V/STOL aircraft 
is dependent on factors other than angle-of-attack stability. For this particular wing 
incidence, the vehicle was statically stable with respect to angle of attack (see fig. 6(a)) 
but the tail-incidence program led to negative (destabilizing) values of the velocity sta- 
bility derivative M u . An inspection of the coefficient E of the quartic characteristic 
equation in appendix A reveals that both angle-of-attack stability (M w ) and velocity sta- 
bility (M u ) can determine the static stability of the aircraft. (A negative sign of the coef- 
ficient E indicates static instability.) For i w = 10°, the unstable value of velocity sta- 
bility was large enough to cause static instability of the aircraft. These results are 
believed to be further substantiated by the free-flight tests of a tilt-wing model as 
reported in reference 2. The model used for the investigation of reference 2 also had 
excessive values of tail incidence at high transition speeds. In spite of the fact that force 
tests showed the model to have static stability with respect to angle of attack, control- 
fixed motions were observed to be aperiodic divergences. Such motions might, of course, 
be taken to be out -of -trim flight conditions. Although no measurements were made of the 
velocity stability derivative, it is believed that a condition existed similar to that for the 
present configuration — that is, instability of the model due to excessive values of tail 
incidence. 

The root locations of the equations of motion for transition flight (except for the 
nonlinear condition of i w = 25°) are plotted on the complex plane in figure 11. The 
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values of the mode ratios 



and 


given in table IV, together with the data of fig- 


ure 11, reveal that as the wing incidence is progressively reduced from 90°, the complex 
roots describing the unstable oscillation of hovering flight move to the left (less unstable) 
and closer to the origin (increase in period). Then, these complex roots separate into 
two real roots defining one stable and one unstable aperiodic mode. The negative real 
root of hovering flight (i w = 90°) has approximately the same value at the low transition 
speeds (iw = 65° and 50°) but combines with the other (unstable) real root to form the 
damped short-period oscillation at i w = 10°. 


The roots describing the dynamic stability of the present vehicle are seen to follow 
a distinct path as the transition to forward flight progresses (fig. 11). The real aperiodic 
modes of hovering flight become the conventional-aircraft short-period oscillation of 
forward flight whereas the unstable oscillation of hovering flight becomes the phugoid 
oscillation of forward flight; although, for this configuration, the phugoid oscillation is 
broken down into two real roots because of static instability with respect to speed. Addi- 
tional calculations for i w = 10° with a value of M u = 0 show that a phugoid oscillation 
results with little change in the short-period mode. In a subsequent section of this paper, 
the path of the roots is shown to depend heavily on the relative values of M w and M u . 
For example, the hovering flight oscillation became the short-period oscillation in for- 
ward flight during flight tests of an experimental tilt -wing aircraft (see ref. 13). 


Correlation of Calculated and Experimental Results 

As previously stated, the model used in the measurement of the stability derivatives 
had undergone a series of free-flight tests during which measurements were made of 
observed control-fixed longitudinal dynamic stability characteristics. A complete 
description of the free-flight technique used for the tests is given in reference 3. A 
photograph of the 1/9-scale model in free flight in the Langley full-scale tunnel is pre- 
sented as figure 12. The dynamic-stability measurements consisted of time histories 
obtained from motion-picture records of control-fixed motions occurring after random 
atmospheric disturbances from trim flight conditions. No attempt was made to match 
calculated time histories because of the random nature of the disturbance input; instead, 
the basic character of the motions in terms of period and damping are compared. Motions 
measured using the free-flight model technique are those which would be seen by an 
observer flying alongside the model at the steady trim speed U Q . Typical control-fixed 
longitudinal motions of this type, based on scaled-up model values obtained during the 
tests, are presented for several wing incidence angles in figure 13. The model motions 
are seen to agree with the results of the calculations (table IV) in an overall sense in that 
the model had an unstable oscillation in hovering flight and became less unstable at high 
transition speeds. The calculated period and damping of the oscillatory mode of hovering 
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flight (P = 9.68 sec, t2 = 2.38 sec) agree well with the scaled-up model values 

(P = 10.4 sec, t2 = 2.5 sec). As the transition to forward flight progressed and the wing 
incidence was reduced to 65°, the oscillation became longer in period (P = 12.4 sec for 
scaled-up model value and 14.11 sec for calculated value) than that for hovering flight. 
Because of the unstable oscillation, the aperiodic divergence could not be confirmed, but 
the model certainly seemed to be simply diverging in a nose-down sense at the end of the 
record. As wing incidence was further reduced to 25°, an unusually long period "oscilla- 
tion" was noted for an initially nose-down motion. This flight condition emphasizes the 
fact that although the time histories of figure 13 are informative, any inspection of limited 
samples of the model motions under certain circumstances may lead to erroneous con- 
clusions regarding the dynamic stability of free-flight models. This is particularly true 
for conditions involving very nonlinear and unsymmetrical aerodynamic characteristics. 
The previously discussed aerodynamic nonlinearity for i w = 25° was not known at the 
time the free-flight tests were conducted; as a result, the time history is based on nose- 
down motion. Additional examination of motion pictures taken during the flights revealed 
that the model would not oscillate when disturbed in a nose-up sense but would pitch up, 
as indicated in figure 10. These motions were regarded as a result of out-of-trim condi- 
tions during the free-flight tests. At the lowest wing incidence of the investigation (10°), 
the model motions observed during the free-flight tests were a short-period oscillation 
in angle of attack and pitch angle and a lightly damped phugoid oscillation or aperiodic 
divergence in displacement. 

The fact that the results of the calculations are in fairly good agreement with the 
motions observed during the free-flight tests indicates that linearized stability calcula- 
tions can be applied to tilt-wing aircraft for most flight conditions. More sophisticated 
time-consuming means of analysis may be required for certain flight conditions in which 
nonlinear and unsymmetrical aerodynamic characteristics prevail. 

Effect of Individual Stability Derivatives on Dynamic Stability 

The values of the various longitudinal static and dynamic stability derivatives were 
varied with the use of the root-locus method in order to determine the effect of the indi- 
vidual stability derivatives on dynamic stability. Calculations similar to those previously 
discussed were also made in which the basic values of the stability derivatives were 
doubled and halved in order to illustrate the magnitude of each derivative at various 
points along the locus. The derivatives were varied at all wing incidences for which the 
dynamic stability calculations previously discussed were made except i w = 25° (the 
nonlinear case). Since preliminary calculations indicated that the dynamic derivatives 
Xq, X w , and Zq were of negligible importance throughout the range of flight conditions 
investigated, the effects of these derivatives were not studied in detail. Consequently, 
results are not presented for these particular derivatives. The following table presents 
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an index to the figures in which the root-locus results are presented for the various 
derivatives: 


Derivative 

Figure 

Hovering flight, i w = 90° 


M u 


14 


X u 


15 


Mq and Mg .... 


16 



Transition flight 



o 

m 

CO 

11 

Jf 

i w = 50° 

t-H 

II 

Jf 

X u 

17 

23 

29 

Zu 

18 

24 

30 

Zw 

19 

25 

31 

M u 

20 

26 

32 

M w 

21 

27 

33 

Mq and . . . . 

22 

28 

34 


Hovering flight.- The locations of the roots presented in figure 14 for the various 
values of M u indicate that the roots of the characteristic equation for hovering flight 
are relatively sensitive to changes in velocity stability. Increases in M u lead to a 
more unstable oscillation with a corresponding decrease in period. It is interesting to 
note that the oscillation becomes more unstable in terms of time required to double ampli- 
tude (that is, a becomes more positive) but the instability remains about constant in 
terms of number of cycles required to double amplitude (o/oj ~ Constant). The locations 
of the roots presented in figure 15 show that although negative increases in Xu tend to 
make the modes of motion less unstable, the roots are relatively insensitive to changes in 
this derivative. Figure 16 shows the location of the roots representing the oscillation for 
pitch-rate stabilization (Mg = 0, Mq = -2 per rad-sec), pitch-attitude stabilization 
(Mg = -1 per rad-sec2, Mq = 0), and various ratios of pitch-attitude stabilization to pitch- 
rate stabilization. With the addition of pitch-rate stabilization only (Mg/Mq = 0), the 
unstable oscillation can be made stable, but the vehicle response to control inputs will be 
sluggish. For the case of artificial stabilization in attitude alone (Mg/Mq = °°), the fre- 
quency of the oscillation is increased with little change in the damping characteristics. 

The most effective means of stabilizing the unstable oscillation is through use of artificial 
stabilization in the form of a combination of pitch-rate (q) stabilization and pitch -attitude 
(6) stabilization. 
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The unstable oscillation of hovering flight has not been of concern to pilots of 
V/STOL aircraft when flying in weather conditions where visual motion cues are avail- 
able inasmuch as the pilot will add the rate and attitude stabilization necessary to sta- 
bilize the aircraft. Some difficulties, however, may arise during hovering flight when 
using instrument references only. Significant improvements in handling qualities were 
obtained by use of artificial stabilization during the flight tests of reference 7. In any 
event, at least artificial pitch-rate stabilization would probably be required for satisfac- 
tory handling qualities for this particular configuration. 

Transition flight .- The root-locus plots presented in figures 17 to 22 for i w = 65° 
show that the stability of the oscillatory mode is very sensitive to changes in M u , Mq, 
and Mg and relatively insensitive to changes in X u and Z u . The unstable oscillation 
could be made stable by the addition of artificial stabilization in Mq or Mg. The 
unstable aperiodic mode is most sensitive to changes in the stability derivatives Z w and 
M w because of the relatively large angle -of-attack content in that mode for this low- 
speed transition flight region (see table IV). 

The data presented in figures 23 to 28 for i w = 50° indicate that the aircraft 
dynamic stability characteristics are sensitive to changes in all the stability derivatives 
except X u and Z w . The unstable oscillation is again stabilized by the addition of 
pitch-attitude or pitch-rate stabilization. The only stability derivative for which changes 
can stabilize the aperiodic divergence without destabilizing the oscillatory mode is M w 
(fig. 27). 

For i w = 10°, the root-locus plots shown in figures 29 to 34 indicate that the sta- 
bility derivatives are beginning to assume the relative importance normally associated 
with conventional aircraft stability derivatives. For example, negative increases in Z w 
tend to add damping to the short-period mode (fig. 31), negative increases in M w tend 
to increase the frequency of the short-period oscillation (fig. 33), and negative increases 
in damping in pitch (Mq) increase the damping of the short-period mode with little effect 
on the phugoid roots (fig. 34) . The unstable aperiodic mode of the basic configuration 
may be stabilized by a negative increase in Z u (fig. 30), a reduction in the negative 
value of Z w (fig. 31), a reduction in the negative value of M u (fig. 32), or a negative 
increase in M w (fig. 33). 


SUMMARY OF RESULTS 

The results of an analytical investigation of the dynamic longitudinal stability of a 
tilt-wing V/STOL aircraft may be summarized as follows: 

1. The results of the calculations using the classical linearized equations of motion 
agree with the results of free-flight tests of a 1/9-scale model of the vehicle. However, 
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the linearized equations may not be applicable to the tilt -wing aircraft for some flight 
conditions in which aerodynamic nonlinearities are present. 

2. The control-fixed longitudinal motions of the tilt -wing aircraft without artificial 
stabilization in hovering flight were dominated by an unstable oscillation similar to that 
displayed by most helicopters. As the transition to conventional forward flight pro- 
gressed, stability characteristics were encountered in which aperiodic divergent modes 
of motion, as well as unstable oscillations, were present. The conventional-aircraft 
short-period and phugoid oscillations began to appear at the high-speed end of the 
transition. 

3. The unstable oscillation occurring in the hovering and low-speed flight regions 
can be stabilized by the addition of a combination of pitch-rate and pitch-attitude stabili- 
zation, but angle -of -attack stability must be increased if the aperiodic divergences are 
to be made stable. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 23, 1967, 

721-01-00-26-23. 
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APPENDIX A 


EQUATIONS OF MOTION 

The linearized, small-perturbation equations of motion for fuselage -level horizontal 
flight (ref. 11), referred to a body system of axes (fig. 1), may be expressed as follows: 

Vertical force: 

-Z u u + (s - Z w )w - (U 0 + Zq)s0 = 0 (Al) 

Longitudinal force: 

(s - X u )u - X w w + (g - Xqs)6> = 0 (A2) 

Pitching moment: 

-M u u - (M w + M^s)w + (s2 - Mqs - M e )6 = 0 (A3) 

For nontrivial solutions, s must be a root of the characteristic equation 

As^ + Bs 3 + Cs2 + Ds + E = 0 (A4) 

where 
A = 1 

B = Mw(-Zq - Uo) - Mq - Xu - Z w 

C = Mw(U 0 X u + Z q X u - XqZ u ) - M W (U 0 + Z q ) + X U (Z W + Mq) - XqM u - Mg + MqZ w - Z U X W 

D = M^(gZu) + Xq(M u Z w - Z U M W ) + Z q (X u M w - M U X W ) + M 0 (Z W + Xu) + M w (U 0 Xu) 

+ X w (MqZ u - UoMu) + gMu - MqZ-yyXu 

E = g(Z u M w - M U Z W ) + Mg(Z u X\y - Z^yX u ) 

The damping and period of a mode of motion, in seconds, are given by the equa- 
tions ti/2 = - Q and P = respectively, where a and u> are the real and 

imaginary parts of the root of the characteristic equation. Additional stability character- 
istics may be obtained by the following relations: 

c l/2 = tl/^/P w n = + “> 2 ? = 5“ 


17 



APPENDIX B 


SIMPLIFIED EQUATIONS OF MOTION FOR HOVERING FLIGHT 

Considerable simplification of the general equations of motion given in appendix A 
can be made for hovering flight. Inasmuch as the stability derivatives X w , M w , and 
Mw will usually be negligible, the vertical-force equation becomes uncoupled from the 
horizontal-force and pitching-moment equations. The resulting characteristic equation 
is 

(s - Z w )js3 - (X u + Mq)s2 + (X u Mq - M u Xq - Mq)s + (X u Mfl + M u g)J = 0 (Bl) 

As a consequence of the uncoupling of the vertical degree of freedom, one root (equal in 
value to Z w ) is immediately known. The remaining cubic describes the oscillation 
involving the fore-and-aft and pitch degrees of freedom. 
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TABLE I.- MASS AND GEOMETRIC CHARACTERISTICS 
OF THE V/STOL TRANSPORT 


Gross weight, lb (N) 

Moment of inertia in pitch, slug-ft^ (kg-m2) 

Wing loading, lb/ft 2 (N/m2) 

Fuselage: 

Length, ft (m) . . 

Cross-sectional area, maximum, ft 2 (m2) 

Height, ft (m) 

Width, ft (m) 

Wing: 

Area, ft^ (m^) 

Span, ft (m) 

Aspect ratio 

Mean aerodynamic chord, ft (m) .... 

Airfoil section 

Tip chord, ft (m) 

Root chord, ft (m) 

Taper ratio 

Sweepback of quarter chord, deg .... 

Dihedral angle, deg 

Pivot location, percent root chord . . . 
Aileron, each: 

Chord, percent wing chord 

Area, ft^ (m2) 

Flap, each: 

Type 

Chord, percent wing chord 

Span 

Slat, each: 

Inboard, 0.45 wing semispan to 

0.69 wing semispan 

Outboard, 0.85 wing semispan to 

1.00 wing semispan 


. 52 000 (231 307) 

. 205 000 (277 940) 

. . . . 97.2 (4 654) 

50.0 (15.2) 

81.8 (7.6) 

12.2 (3.7) 

9.1 (2.8) 

. . . . 534.4 (49.6) 

. . . . 67.5 (20.6) 

8.53 

8.07 (2.5) 

. . . NACA 633-318 

6.0 (1.8) 

9.8 (2.9) 

0.61 

4.13 

-2.12 

42.5 

25 

31.1 (2.9) 


Double slotted 
...... 47 

..... Full 


Chord, 0.20 wing chord inboard 
to 0. 10 wing chord outboard 

. . . Chord, 0.10 wing chord 
full length 
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TABLE I.- MASS AND GEOMETRIC CHARACTERISTICS 
OF THE V/STOL TRANSPORT - Concluded 


Vertical tail: 

Area, ft 2 (m 2 ) 130.0 (12.1) 

Span, ft (m) 15.6 (4.8) 

Aspect ratio 1.87 

Airfoil section: 

Root NACA 0018 

Tip NACA 0012 

Tip chord, ft (m) 3.3 (1.0) 

Root chord, ft (m) 13.3 (4.1) 

Taper ratio 0.25 

Sweepback of quarter chord, deg 26 

Rudder: 

Tip chord, ft (m) 1.1 (0.3) 

Root chord, ft (m) 3.8 (1.2) 

Tail length, center of gravity to 0.25 mean aerodynamic chord, ft (m) . . . 21.4 (6.5) 

Horizontal tail: 

Area, ft 2 (m2) 170.9 (15.9) 

Aspect ratio 5.68 

Airfoil section: 

Root NACA 0015 

Tip NACA 0012 

Tip chord, ft (m) 3.5 (1.1) 

Root chord, ft (m) 7.0 (2.1) 

Span, ft (m) 31.14 (9.5) 

Taper ratio 0.50 

Sweepback of quarter chord, deg 9.50 

Mean aerodynamic chord, ft (m) 5.5 (1.7) 


Tail length, center of gravity to 0.25 mean aerodynamic chord, ft (m) . . . 24.8 (7.6) 


Propellers: 

Main: 

Number of blades 4 

Diameter, ft (m) 15.5 (4.7) 

Tail: 

Number of blades 3 

Diameter, ft (m) 8.0 (2.4) 

Moment arm, wing pivot to rotor center, ft (m) 32.0 (9.8) 
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TABLE H.- SCALING FACTORS 


To scale model values up to full-scale values, X 


l 


Dimensional and mass properties 


To scale 

Linear acceleration . 

Length 

Area 

Mass 

Force 

Moment 

Moment of inertia 
Linear velocity . . . 

Time 

Angular velocity . . 
Angular acceleration 


Stability derivatives 


To scale 


Multiply by 

1 

X 

X2 

X3 

X3 

X4 

X5 

x !/2 

xl/2 

X-l/2 

,-l 


Multiply by 





TABLE EL- STABILITY DERIVATIVES 


X W) 

Xq, 

per sec 

ft/sec-rad 


(a) U.S. Customary Units 


2ii i 


90 

0 

-0.2855 

0 

65 

43.5 

-.1848 

0 

50 

51.9 

-.1218 


25 

(Nose up) 

108i 3 

-.1058 


25 

(Nose down) 

108.3 

-.1058 



(b) SI Units 


deg m/sec per sec per sec m/sec-rad per sec | per sec | m/sec-rad | per m-sec | per m-sec | per rad-sec 




0 

-.1940 

-.2235 

-.3553 

-.3553 


90 

0 

-0.2855 

0 

65 

13.3 

-.1848 

0 

50 

15.8 

-.1218 


25 

(Nose up) 

33.0 

-.1058 


25 

(Nose down) 

33.0 

-.1058 



.0089 


-.0105 


-.7213 





















TABLE IV.- SUMMARY OF RESULTS FOR BASIC AIRCRAFT 


V’ 

deg 

90 

Mode 

Oscillatory 

Aperiodic 

Aperiodic 

O’ 

0.29 

0 

-.87 

w 

±0.65 

0 

0 

*1/2’ 

sec 

(*) 

-2.38 

OO 

.80 

P, 

sec 

9.68 

c l/2 

(*) 

-0.246 

w n 

0.711 

e 

(*) 


w 

u 



e 

u 


deg 

deg 




ft/sec 

m/sec 

-0.410 

0 

OO 

0 

1.55 

1.04 

5.09 

3.41 

65 

Oscillatory 

0.16 

±0.44 

-4.29 

14.11 

-0.304 

0.474 

-0.341 

0.59 

0.99 

3.25 


Aperiodic 

.08 

0 

-8.99 





2.33 

.46 

1.51 


Aperiodic 

-.78 

0 

.89 






.96 

1.09 

3.58 

50 

Oscillatory 

0.04 

±0.32 

-16.50 

19.39 

-0.851 

0.327 

-0.129 

0.41 

0.69 

2.26 


Aperiodic 

.39 

0 

-1.77 





1.33 

1.21 

3.97 


Aperiodic 

-.86 

0 

.81 





1.88 

1.92 

6.30 

25 

Oscillatory 

-0.09 

±0.24 

8.09 

26.67 

0.303 

0.251 

0.342 

0.41 

0.38 

1.25 

(Nose up) 

Aperiodic 

.48 

0 

-1.45 





3.43 

1.99 

6.53 

. 1 

Aperiodic 

-.84 

0 

.83 





5.81 

2.87 

9.42 

25 

Oscillatory 

-0.28 

±0.36 

2.48 

17.57 

0.141 

0.454 

0.616 

2.52 

1.33 

4.36 

(Nose down) 

Oscillatory 

.01 

±.14 

-46.70 

44.91 

-1.040 

.141 

-.105 

1.40 

.59 

1.94 

10 

Oscillatory 

-0.60 

±0.75 

1.15 

8.42 

0.137 

0.959 

0.628 

7.29 

2.05 

6.73 


Aperiodic 

.09 

0 

-7.96 






.68 

.42 

1.38 


Aperiodic 

-.15 

0 

4.74 






.86 

^02j 

.07 


^Negative signs indicate unstable modes of motion. For example, if 
1 1/2 = -2.38, then t2 = 2.38. 
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Figure 5.- Variation of forces 
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Figure 9.- Location of the roots on the complex plane for hovering flight. 
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Figure 12.- Photograph of the 1'4-scale model in free flight. 
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Figure 19.- Root locus for Z w . i w = 65°. 
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Figure 22.- Root locus for M a and Mq. i w = 65° 
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Figure 26.- Root locus for M| 
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Figure 33.- Root locus for M, 
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